Elliptic billiard 



— a non-trivial integrable system 
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We investigate tiie semiclassical energy spectrum of quantum elliptic billiard. The nearest neigh- 
bor spacing distribution, level number variance and spectral rigidity support the notion that the 
elliptic billiard is a generic integrable system. However, second order statistics exhibit a novel prop- 
erty of long-range oscillations. Classical simulation shows that all the periodic orbits except two 
are not isolated. In Fourier analysis of the spectrum, all the peaks correspond to periodic orbits. 
The two isolated periodic orbits have small contribution to the fluctuation of level density, while 
non-isolated periodic orbits have the main contribution. The heights of the majority of the peaks 
match our semiclassical theory except for type-O periodic orbits. Elliptic billiard is a nontrivial 
integrable system that will enrich our understanding of integrable systems. 



I. LEVEL STATISTICS OF ELLIPTIC BILLIARD 

The boundary of elliptic billiard (EB) is defined as 

la^ + 2/^/^^ = 1 with aspect ratio a = b/a and ah — 1 
There arc four symmetry classes of eigenvalues 
and eigenfunctions according to the reflection symmetry 
of eigenfunctions against the x and y axes. We study 
the level statistics of the odd-odd class, which implies 
a quarter EB with Dirichlet boundary condition. The 
eigenvalues are computed by direct diagonalization of the 
Hamiltonian [ll-Q- 

There are two kinds of averaging used for analysis of 
statistics in quantum chaos: ensemble averaging, e.g., 
over realization of disorder, and spectral averaging Q. 
For integrable systems, only spectral averaging had been 
done before Refs. @, 0], in which the system parameter 
is sampled from algebraic numbers around the central 
value. A better way of sampling is sampling from a nor- 
mal distribution |8|. In this article, for EB a is sampled 
from a normal distribution centered around 1/2 and for 
circular billiard (CB) around 1. 

The nearest neighbor spacing distribution P{s) with 
the level spacing s, level number variance I](e,i?) with 
energy e and energy interval width E, and spectral rigid- 
ity /S.z{e,E) are shown in Fig. [TJ P{s) generally follows 
Poisson distribution, but displays deviation at small s. 
This is a general phenomenon showing level repulsion of 
integrable systems Q. The oscillations of Yi{e,E) are 
explained by the semiclassical equation 
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where j counts periodic orbits (PO), h is the Planck con- 
stant, Aj^ Tj amplitude and period of PO respectively, 
2N the dimension of phase space Like rectangular 
billiard (RB), the match between numerical and theo- 
retical results is achieved only after ensemble averaging 
in theoretical calculation Spectral rigidity saturates. 
Like RB, A3 of EB generally follows the squared-energy 
scaling relation Q : 



but displays long-range oscillations. All the statistics 
support that EB is a generic integrable system. 



A. spectral rigidity and global variance 

The long-range oscillations of spectral rigidity show 
some global oscillations of level density. To reveal this 
global oscillations, we investigate a special case of level 
number variance, the global variance defined as 

^,ie) ^ {[6^{e)Y) ^ ([^(6) - (^(6)>]^>, (3) 

where ^{e) is the spectral staircase. Global variance is 
also a special case of the correlation function of spectral 
staircase defined as {5.yV {ei)6.yV {€2)) It is proved 
that the the rigidity and global variance are equal: 0, @ 
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The above equation is partially correct as Sg fluctuates 
around A3 for RB, CB and EB as demonstrated in Fig. 
[21 For RB, Sg and A3 follow the squared energy scaling 
relation. For CB or EB, the long-range oscillations of Sg 
and A3 are "synchronized" . The scales of oscillations of 
S(e, E) and A|^(e) are different. The former is i? ^ -y/e, 
and the latter is longer than this. 

We think the long-range oscillations of A3 and Eg orig- 
inate from collective effect of type-R POs of EB or CB 
and exist in any billiard with whispering gallery modes. 
The whispering gallery POs have approximately equal or- 
bit length or multiples of that and hence they coherently 
cause global fluctuation of level density. 



II. PERIODIC ORBITS OF ELLIPTIC 
BILLIARD 

Semiclassical spectral properties are closely related to 
classical POs. We have the following theorems concern- 
ing classical POs of EB, demonstrated in Fig. [3l First, 
the numerical work of hundreds of families of POs veri- 
fies that all the POs except two on axes are continuous. 
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FIG. 2: Comparison between spectral rigidity and global vari- 
ance. Red line: A3°(e). Blue line: 'Sg{e) of RB, quarter EB, 
and quarter CB. The ensemble size of RB calculating Eg is 
10^. The oscillation of Sg(e) is unlikely to disappear in the 
limit of infinite ensemble size. 



FIG. 1: EB's P{s) (a), E (b), A3(e,-E) saturation with E (c), 
long-range oscillations of A^(e) with energy interval width 
(d), calculated by averaging over an ensemble of 1350 sam- 
ples of a. (a) Black dots: numerical result. Red line: Poisson 
distribution. Inset: shape of quarter EB. (b) Black line: nu- 
merical result. Green line: theoretical result calculated by 
averaging over variance calculated from 200 shortest POs. 
(c) Black line: numerical result. Red line: saturated value 
of As{e,E). (d) Black line: saturated A3°(e). Red line: the 
scaling relation in Eq. (2] 




Inside the area covered by a family of POs, there are 
four directions to find POs in the family. This makes 
the factor c in Eq. [12] below same for every continu- 
ous families of POs. Second, the envelope of a fam- 
ily of POs is an ellipse or a hyperbola. Third, the 
envelope ellipses, hyperbolas and EB arc confocal. In 




FIG. 3: Rs,!, Ra,i, Rb,2, and O4 families of POs. 
envelope ellipse or hyperbola. 



Red line: 



Fig. [3l the foci of the envefope eUipse for R31 are at 
X = ±Vl. 2282682 - 0.092972 = ±1/372, which is same 
as EB; the foci of the envefope hyperbola of O4 are at 
X = ±^1.15472 + 0.4082482 = ±vW2- Fourth, all the 
POs of a family with the same envelope curve have the 
same orbit length, for example all the Rs^i POs have or- 
bit length 6.0322. Fifth, for type-R POs, every orbit line 
is tangent with the envelope ellipse and for type-0 ei- 
ther the orbit line or its extended line is tangent with 
the envelope hyperbola. 

III. FOURIER ANALYSIS OF SPECTRUM OF 
RECTANGULAR, CIRCULAR BILLIARDS 

To reveal the underlying POs of semiclassical spec- 
trum, we perform Fourier analysis on the spectrum. In 
this section and the section below, we do not perform 
ensemble averaging on the Fourier analysis. The posi- 
tions of peaks of Fourier analysis are identified with or- 
bit lengths. Numerical work of EB confirms that all the 
peaks have the same half-width. Hence the height of a 
peak is proportional to amplitude of a peak and propor- 
tional to the amplitude of a PO. 

height of peak on amplitude of peak cx amplitude of PO. 

(5) 

A. Fourier analysis of spectrum 

The Fourier analysis of spectrum has been success- 
fully applied to analyze Sinai billiard, hydrogen atom 
in a strong magnetic field [l^l- The trace formula ex- 
presses the fluctuation of level density as a summation 
over (families of) POs: 

p(fc)=^^,(L,)e'^^^+", (6) 
j 

where Lj is orbit length of a family of POs and a the 
Maslov index. We do not consider a in this article. If we 
extend Aj{Lj) to A{1): the amplitude of POs over orbit 
length by extending Lj to the whole line of I, we have 

/oo 
A(/)e'^'dL (7) 
-00 

The reverse of a PO has orbit length —I. The inverse 
Fourier analysis gives the amplitude 

A{1) = —j p{k)e-^'^dk. (8) 

The discreteness of spectrum makes Fourier analysis into 
a summation over eigen- momentum fc^, where i counts 
eigenvalues: 

^(0 = iE^"''-'' (9) 
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For simplicity, we ignore the subscript j below. 

B. rectangular billiard 
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FIG. 4: Black line: Fourier analysis of RB with two sides as 1 
and (^/5^- l)/2. Red square: semiclassical theory. The L — 2 
peak is normalized to 1 in both numerical and theoretical 
calculations. 

A direct quantum mechanical calculation proves that 
for RB the amplitude of a family of POs is given by 

A(L)cx-^, (11) 

where the momentum space factor c = 1/2 for POs par- 
allel to boundary, and 1 otherwise @, Q as for the former 
there are two directions to find POs from the same family 
while for the latter there are four directions. Although 
Eq. [TT] is simple, it gives correct position and height 
of almost every peak shown in Fig. |4l The rare large 
discrepancy for some POs is caused by interference of 
close-by POs. 

C. circular billiard 

For CB, different POs cover different areas in position 
space. To account for this difference, we add the factor 
5', the area covered by a family of POs in position space 
(in the case of EB, the area covered by orbit lines in Fig. 
[J) and the amplitude of POs reads 

rS* 

A{L) (X (12) 

where the momentum space factor c = 1/2 for POs along 
diameters, and 1 otherwise as for the former there are two 
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Circular billiard 




FIG. 5: Black line: Fourier analysis of CB. Red square: semi- 
classical theory. The Z/ = 4 peak is normalized to 1 in both 
numerical and theoretical calculations. 



directions to find POs from the same family while for 
the latter there are four directions. For a family of POs 
denoted by (n,m) with m < n/2, where in a period, the 
particle collides n times with the boundary and makes m 
revolutions, 



Given its similarity with Gutzwiller trace formula [llj . 
this formula gives some new information compared with 
Berry- Tabor trace formula (T^j. It is worthwhile to un- 
derstand its relation with the latter. 



IV. 



FOURIER ANALYSIS OF SPECTRUM OF 
ELLIPTIC BILLIARD 



Elliptic billiard 




L ^ 2nsin(^). 



(13) 
(14) 



Eq. [12] is verified in Fig. [5] except for interfering close- by 
orbits. 



FIG. 6: Red arrows: positions of unstable and stable isolated 
POs. The Fourier analysis contains eigenvalues of all the four 
symmetry classes of EB with a — 1/2 but the almost degener- 
ate eigenvalues are removed by only considering one of them. 
The heights of peaks are normalized by defining the peak for 
iig.i equal to 1 in both numerical result and theory. 



D. 



summary 



The position and height of peaks of RB and CB are 
correctly given by the amplitude formulae of POs in Eqs. 
[m and [121 From our knowledge, the factor S in Eq. [T2| 
has never been pointed out before. Its existence in EB 
is also confirmed below. An argument to justify Eqs. 1111 
and [T2l is as follows. 

The amplitude of a family of POs is the linear super- 
position of all the individual POs in the family. For inte- 
grable systems, we assume that the amplitude of each PO 
is only decided by its orbit length and given by 1/ \/Z, 
which is derived theoretically in the case of RB or more 
generally from Berry- Tabor formula . In this assump- 
tion, orbit stability plays no role. All the POs of a fam- 
ily contribute equally to the amplitude. The factor cS is 
used to quantify how many POs are in a family, c gives 
the number of POs in the momentum space and S in the 
position space. For RB, every family of POs covers the 
whole billiard and S plays no role therein. 

Now we have a trace formula for integrable systems 
from Eq. [6) 
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Comparison of Fourier analysis of EB with the semi- 
classical theory in Eq. [T^] is shown in Fig. [HI The po- 
sitions of peaks match well with very small discrepancy. 
At L = 2.83, the short peak is due to the unstable PO 
with L = 46, indicating that the unstable isolated PO is 
unimportant to the level fluctuation. At L = 5.66, the 
stable isolated PO with L = ia contributes another short 
peak. The isolated POs contribute little to level fluctu- 
ation, which was already demonstrated in the modified 
Kepler problem containing one stable isolated PO Q. 
The unimportance of isolated POs can be explained by 
our argument in the last section as its position space 
factor S = 0. At L = 6.32, the continuous Ri^i-POs 
contribute significantly. Around L = 6.6, there are many 
peaks overlapping together. These are the contribution 
from Rn.i with n ^ 1. These POs has the orbit length 
around the perimeter cc^ual to 6.85. 

For heights of peaks, good match between numerical 
and theoretical results is achieved for type-R POs except 
for i?5^2 but not for type-0 POs. Two kinds of discrepan- 
cies exist. First, for O4, Og, Oio, O14, Oig, the discrep- 
ancy is a factor around 1. Second, for Og the discrepancy 
is very large. For the equally spaced 04,204,304,404 
(repetition of O4) and 06,206,306 (repetition of Oq), 
their relative magnitudes are correctly given by the fac- 
tor l/VX. 
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V. CONCLUSIONS 

Although EB is a generic integrable system, it contains 
several new characteristics. The second order statistics 
including spectral rigidity and global variance display 
long-range oscillations, not noticed in other systems be- 
fore. We think such oscillations exist for any billiard with 
whispering gallery modes. The global variance seems to 
be a statistical instrument worth further studying. The 
Fourier analysis of spectrum gives the position and am- 



plitude of POs. For RB and CB, the numerical results 
match our proposed theory of amplitude of POs. For 
EB, only the type-R POs fully match our theory. Fur- 
ther work is needed on type-0 POs. Future work will ad- 
dress the orbital magnetic response of EB. Typc-R and 
type-0 POs enclose different magnetic flux and hence 
they should show different orbital magnetic response. In 
summary, EB is a non-trivial integrable system that will 
further enrich our understanding of quantum chaos. 
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